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Abstract: Einstein-Maxwell theory implies the mixing of photons with gravitons in an 
external electromagnetic field. This process and its possible observable consequences 
have been studied at tree level for many years. We use the worldline formalism for 
obtaining an exact integral representation for the one-loop corrections to this amplitude 
due to scalars and fermions. We study the structure of this amplitude, and obtain exact 
expressions for various limiting cases. 



1. Photon-graviton mixing at tree level and one loop 

As has been realized already in the sixties [1] Einstein-Maxwell theory in a constant 
electromagnetic field contains a tree level vertex for photon-graviton conversion, 



KhrFT* - ^F^fafi- (1) 

Here h^ u denotes the graviton, f^ v the photon, and F^ u the external field, k is the 
gravitational coupling constant. This vertex implies the possibility of photon-graviton 
oscillations [H El El IH 0, EJ [7] which are analogous to the better-known photon-axion 
oscillations in a field [5]. In the presence of an external field, the true eigenstates of 
propagation can be obtained by solving the following system of dispersion relations [3] 



(rpPtf - k a k? \kC kX ^ \ ( a p (k) \ _ 

\ freC^ ^(r]^ K r] vX +r]^ x rj UK -2rj^T] KX +...)) \h KX (k) ) ~ ' ^ 

Here C^ u ' a denotes the Fourier transform of the vertex Q, 



C u ' a = (F ■ fc)V" + F" a k u + F ua k» - (F ■ kfr] va - (F ■ k) u r]^ a . (3) 

For small deviations from the vacuum dispersion relations k 2 = this second order 
equation can be linearized. An efficient formalism for solving the system ^ under this 
condition was developed in [3]. 

Taking one-loop corrections into account, the dispersion relation matrix gets 
modified in the following way (9j ITU] , 



fr} afi k 2 - k a k? - Ii a ^ {kC kX '° - U^ a \ ( a p (k) \ _ 

V \kC^ -W- v ^ ^(7]^ K 7] uX +r]^ x r] UK -27]^ u r] KX +..)-YI^' kX )\h KX (k) J ~ ' 

(4) 

Here Yl a ^, IP 1 ''' 3 , and fl flu,KX denote the one-loop photon-photon, graviton-photon, and 
graviton-graviton vacuum polarization tensors in a constant field. These quantities in 
principle have to be calculated with all possible loop particles. Eqs. Q generalize the 
QED dispersion relation 



(r] a ^k 2 - k a k? - IF ,/3 ) ap(k) = . (5) 

This case is well-known and has been studied by many authors (see, e.g., [HI H21 [13] ) . It 
leads to a complicated dependence of the phase velocity on polarization, field strength 
and frequency (see [14J for a detailed discussion). 
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In this talk, we report on our recently concluded calculation [91 [10] of the photon- 
graviton polarization tensor in a constant electromagnetic field II' u/ ' a , with a charged 
scalar or spin | particle in the loop. As a Feynman diagram, this amplitude is 
represented by fig. [Tj 



h 



-4, 



Figure 1. One-loop photon-graviton amplitude in a constant field. The double line 
represents the propagator of a charged scalar or spin ^ particle in a constant field. 

In [9] the worldline formalism was used to obtain compact parameter integral 
representations for this amplitude. The numerical and structural analysis has been 
concluded only recently [10]. Since this formalism is presently still somewhat novel, 
particularly in applications to gravity, we will start with shortly reviewing its basics 
from a user's point of view. 

2. Worldline formalism in a constant electromagnetic field 

The worldline formalism goes back to Feynman's representation of scalar [15] and spinor 
[16J QED in terms of relativistic particle path integrals. Let us write down Feynman's 
integral for the simplest possible case, the one-loop effective action in scalar QED: 



T(A) = ( d*x C(A) = I" ~ e~ m2T [ Vx(r) e 

J JO T Jx(T)=x(0) 



-S[x(t)] 



(6) 



10 T Jx(T)=x(0) 

Here m and T denote the mass and proper time of the loop scalar. The worldline 
path integral is to be calculated over loops in spacetime with fixed periodicity T, and a 
worldline action given by 



S[x(t)] 



dr 



x 2 + ieif AJx{t)) 



(7) 



In the fermion QED case, a number of different ways have been found to implement the 
spin in the worldline path integral. Feynman's original formulas [16] are based on a spin 
factor involving Dirac matrices, but for the purpose of analytic calculation it is usually 
preferable to implement spin by the following additional Grassmann path integral [T7] . 



Vij){f) exp - £ dr ■ i> - ieVF^V 
Here the functions ip^ are anticommuting and antiperiodic, 



(8) 



^(ti)V(t 2 ) = - ^(t 2 )V(ti), VCH = -V(0) 
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(9) 



See, e.g., chapter 3 of [18J for a derivation of the path integral representations ([6]), (j8l) 
from quantum field theory. During the last fifteen years various efficient methods have 
been developed for the evaluation of this type of path integral. We are concerned here 
with the so-called "string-inspired" approach [T9J, [20j |21] (see [18] for a review) which 
aims at an analytic calculation of the worldline path integral (see the contributions by 
G.V. Dunne and K. Klingmiiller to these proceedings for alternative approaches). This is 
achieved by manipulating the path integral into gaussian form, usually by a perturbative 
or higher derivative expansion, and then performing the gaussian integration formally 
using worldline correlators. Those worldline Green's functions are, for the coordinate 
path integral, 



(r -r V 

(^(r 1 K(r 2 ))= -G B {r U T 2 )Tr, G b (t 1 ,t 2 ) = \t 1 ~t 2 \- ' 



T 



(10) 



and for the spin path integral 



(r(riW(r 2 )) = G F {n, r a ) rf\ G F {n, r 2 ) = sign(r 1 - r 2 ). (11) 

We will often abbreviate Gb(ti, t 2 ) =: Gbyi etc. The coordinate Green's function is not 
unique, since it depends on the zero mode fixing of the path integral [18]. The choice 



of ( 10 ) corresponds to a definition of the zero mode as the loop center-of-mass, 



4 := I [ T drx^r). (12) 



o 



T 

To obtain, say, the one-loop N photon amplitude, one expands the Maxwell field in N 
plane waves with given polarization vectors ef 



N 

A»{x{t)) = e?e ik *- x(T) . (13) 



i=i 



This leads to each photon being represented by a photon vertex operator, 

VJL[k, e] = [ T dre- x(r) e lk ^ x{r) (Scalar QED) , 

Jo 

V S p in [k, e}= J dr[e- x(t) + 2ie ■ipk-tp] e ik ^ xir) (Spinor QED) . 

(14) 

After expanding the interaction exponential to iVth order, the path integrals are gaussian 



and can be evaluated by the correlators ( 10 ) , (11). This leaves one with the global proper 



time integral, and one parameter integral for each photon leg. 
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It has emerged that this formalism is particularly well-suited to the calculation of 
QED amplitudes in a constant background field [221 123] • The reason is that, once one 
has obtained the parameter integral representation for a given amplitude for the vacuum 
case, one can construct the corresponding integrals in the presence of a background field 
with constant field strength tensor F^ v by the follwing simple substitutions 

• Change the worldline Green's functions: 



G B (n,T 2 ) -> Gb(t 1 ,t 2 ) 



-^— 2 ( f— e -fag"k» + ieFG Bl2 - ^ 



-ieFTG B12 



cos(eFT) 



(15) 



(A 'dot' on a Green's function denotes a derivative with respect to the first variable.) 
Change the free path integral determinants: 



(4ttT)" 



(4ttT)" 



(4ttT)-^ det"5 



(4ttT)' 



D _ 

2 det 



sin(eFT) 
eFT 

tan(eFT) 
eFT 



(Scalar QED) , 



(Spinor QED) . 



(16) 

It should be remarked that the worldline formalism is closely related to the standard 
Fock-Schwinger proper-time representation of propagators in external fields [2H 125] . 
Therefore the resulting integral representations have generally the same structure as 
the ones obtained by that method (see, e.g., [22] and V. Skalozub's contribution to 
these proceedings). However, the worldline approach is more global in the sense that it 
applies directly to a whole loop, rather than to the individual propagators making up 
the loop. This also implies that the worldline integral representations can be written 
down without fixing the ordering of the external legs along the loop. Another advantage 
is that the use of the so-called "Bern-Kosower substitution rule" [20] provides a simple 
way of inferring the spinor loop integrands from the scalar loop ones. This effectively 
circumvents the usual Dirac algebra manipulations. 

In flat space, the gaussian integration of the worldline path integral can be done 
naively, and no ill-defined expressions are produced. For applications to gravity, we need 
to generalize the path integrals (|6]), ^ to gravitational backgrounds. Here we enter 
the realm of path integrals in curved spaces, a subject notorious for its mathematical 
subtleties. 

Naively, one would introduce background gravity by a simple replacement of the 
kinetic term, 
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S = - [ T drx 2 ^\\ T dTx»g, u (x(T))x» . (17) 
4 Jo 4 Jo 

After the usual linearization g^ v = rj^ + kH^ this would lead to a graviton vertex 
operator of the form 

£a V C drx»x v e lk ' x . (18) 
jo 

However, using this vertex operator in a naive gaussian path integration immediately 
leads to ill-defined expressions involving, e.g., <5(0),<5 2 (rj — t/),... A complete 
understanding of these difficulties, and of the steps which have to be taken to solve them 
in the "string-inspired" framework, has been reached only recently [27J Ell ESI G2B E] . 
Here we can only briefly sketch the correct procedure for the spinless case; all the 
necessary details and the generalization to spin half can be found in [32J. 

(i) In curved space, the path integral measure is nontrivial. Exponentiate it as follows, 
Vx = Dx J] v /det^(x(r)) = Dx f DaDbDc e -VM.M ? ( 19 ) 

0<t<T JPBC 

with a ghost action 

S gh [x,a,b,c] = f T dr ^ V (x)(aV + &V). (20) 
Jo 4 



This modifies the naive graviton vertex operator (18) to 



V£jk, e] = e„ v f dr [x»x v + a^a" + Wc v ] e ik ' x . (21) 

J 

(ii) The correlators of these ghost fields just involve 5 functions, 



(b^ny(T 2 )) = - A5(n - T2 )rT ■ (22) 

The ghost field contributions will cancel all divergent or ill-defined terms. 

(iii) This cancellation of infinities leaves finite ambiguities. From the point of view of 
one-dimensional quantum field theory, we are dealing here with an UV divergent but 
super-renormalizable theory which requires only a small number of counterterms 
to remove all divergences. The coefficients of the counterterms have to be fixed in 
a way which reproduces the known spacetime physics. 
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(iv) As it turns out, these counterterms are regularization dependent, and in general 
noncovariant. This points to a violation of covariance by the regularization method. 
Presently, the only known covariance-preserving regularization method is one- 
dimensional dimensional regularization [30l [3Tj , which has only a single covariant 
counterterm proportional to the curvature scalar (— R/4 in the present notations). 

(v) The zero mode fixing leads to further subtleties. The simplest possibility would 
be to fix a point x on the loop, x(r) = x + yij). It leads to the so-called 
DBC ("Dirichlet boundary conditions") propagator for the coordinate field which is 
known to yield the same effective Lagrangian as would be obtained also by using the 
standard heat kernel expansion. For flat space calculations, the "string-inspired" 
choice (12) is generally more convenient since it is the only one which leads to a 
worldline propagator for the coordinate field depending only on t\ — t-i- It can be 
easily shown that the effective Lagrangians obtained in both ways differ only by 
total derivative terms [33J. 

This continues to be true in curved space, however here those total derivative 
terms turn out to be noncovariant in general, with only the DBC choice yielding 
a manifestly covariant form of the Lagrangian. The noncovariance of the total 
derivative terms present in the "string-inspired" approach poses no problem in 
principle but in practice, since it invalidates the use of the Riemann normal 
coordinate expansion, which is an almost indispensable tool for this type of 
calculations. This remaining problem was solved in [31]. There it was shown that, 
using Riemann normal coordinates from the beginning and performing a BRST 
treatment of the symmetry corresponding to a shift of Xo, the difference between 
the two effective Lagrangians can be reduced to manifestly covariant terms. This 
is achieved by the addition of further Fadeev-Popov type terms to the worldline 
Lagrangian in the "string-inspired" scheme. Those terms are infinite in number 
but easy to determine order by order. 

All this can be generalized to the spin half case [32] . Thus a standard scheme of 
calculation is now available for one-loop effective actions and amplitudes involving 
scalar or spinor loop particles and background gravitational fields. See [35] for some 
applications to effective actions and anomalies, |37] to graviton amplitudes. More 
recently also worldline path integrals representing vector and antisymmetric tensor 
particles coupled to background gravity have been constructed [38] . 
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3. Photon-graviton polarization tensor in a constant field 



Returning to the one-loop photon-graviton amplitude in a constant electromagnetic field 
Fy, v , we will now sketch its calculation for the scalar loop case. According to the above, 
this amplitude can be represented by the following expression, 



^n^ai (k)e c 



ien r°°dT 



4 Jo T 



where Z^ = eTF^ u and V^f are the photon and graviton vertex operators (|14h , (1211). 



e- m i (4vrT)-2det-2 



sm(Z) 
Z 



Kscalt^J £ ^]^s4,l[ k,E c 



(23) 



Using the Wick contraction rules (10), (22) yields (Qbu '■= Qb\2 — Qbw etc.) 



n^cai (k) 



seal 



en 



dT 



e- m2T det~5 



sin(Z) 
Z 



T d n f T dr 2 e- k ^ k l^r, 

o Jo 



[ggu - 2 w) (k ■ Gbi2) - \g& 2 ( Gbu ■ ky + (fj, <- v 

+ (Qb\2 ■ k)f{QB\2 ■ k) V {k ■ Q B \2 



(24) 

The T integral has an UV divergence at the lower limit. Using dimensional 
regularization, this divergence takes the form 



scal,div 



ie 2 n 



1 



3(4vr) 2 £> -4 



(25) 



where C' ll/,a is the tree level vertex (J3J). Adding the corresponding counterterm yields 
the renormalized vacuum polarization tensor IT^" obeying the usual renormalization 
condition S££T(Jfe = 0) = 0. 

Next, appropriate photon and graviton polarizations have to be selected, where it 
turns out to be convenient to use the photon polarization vectors also to construct the 
graviton polarization tensors: 

Photon: £_i_,£||, 

Graviton: e®**" = e ±fl e ±u - s^e^, e®"" = e ±fl e^ + e^e ±u . 

Here we have assumed a Lorentz system such that B and E are collinear, and the 
subscripts on the photon polarization vectors refer to the same direction. Further, no 
information is lost by assuming that the photon propagation is perpendicular to the 
field [TO]. 
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With these conventions, the components of the tree level amplitude become 



- 2Buj, 
2Eu, 

- 2Euo, 

- 2Bu. 



(26) 

Here uj = k° = \k\ denotes the photon/graviton energy. Finally, it is convenient to 
normalize the loop amplitude by the tree level one, making the amplitude dimensionless: 



U^(uj,B,E) = Re\ 



n£i(a,i?,s) 
-Uc Aa 



(27) 



(A — ®,®, a =_L, ||). Here we have further introduced the dimensionless variables 



Cj = » B = , E 



eE 

2 



The spinor loop calculation proceeds completely analogously, just with some 
additional terms coming from the evaluation of the spin path integral (|8]). 

At this stage, the four independent components of the scalar or spinor loop 
amplitude are given in terms of two-parameter integrals, with integrands involving 
trigonometric functions of the proper times and external parameters. Let us write 
down here these integrals for the case of a spinor loop and a purely magnetic field: 



n^ n (u,B) = aRej o y e" 



dvfc?(s,v,uj,B) 



(28) 



TT, 



spin 



TT, 



spin 



e|| 

spin 
spin 



r (A 



.M z _ 

47T ltanh(2;) 

{S B12 f - {S F12 f + {A F12 f - (A m2 + A FU ) (A B12 + - + A FU 



+Ab12\(Sb12) 2 — (S F i 2 ) 2 — (Abu + A F n) 2 + (Api 2 ) 2 — vSb12 + S F 12 



2B 



+ 



3 J ' 



An ltanh(; 



exp 



r (A 



B12 



,\ UJ 



2B 



x 



vS B i2 - S F12 ( A B12 + A FU ) + A B \ 2 \vS B \2 - S F12 + l-v 



2 N w_ 

>2B 



+ 



0. 



(29) 
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Here s = —im 2 T, z = iBs, and the integrand involves the standard worldline coefficient 
functions [2"3] 



Abu 



sinh(zt> ) 

sinh(^) ' 
cosh(^v) 1 

sinh(z) z' 



A B \i = A B 22 = coth(z) - 1 , 

-zaa cosh^-u) — cosh(z) 

Ab\2 = A B12 — A Bn = — 



S 



F12 



F12 



sinh(z) 
cosh(zt>) 
cosh(z) 
smh(zv) 



cosh (z) 
A F ii = A F2 2 = tanh(z) . 

(30) 

The parameter v is related to the original proper-time variables r 12 by v — 1 — 2tl/T 
(the translation invariance of the worldline correlators has been used to set r 2 = 0). 
See [10] for the scalar loop and general constant field cases. 

4. Properties, special cases 

Let us now discuss some properties and limiting cases of the amplitude: 

Ward identities: The gauge Ward identity for this amplitude gives the familiar 

transversality in the photon index, 



k a U^ a (k) = 0. (31) 

The gravitational Ward identity, derived from invariance under infinitesimal 
reparametrizations, connects n MI/,a with the corresponding photon-photon polarization 
tensor U^' a (k), 

k,W^(k)= l 2 KF\W> a (k). (32) 

(Similarly, non-transversality was recently found for the gluon polarization tensor in a 
chromomagnetic background field [39] .) 

Selection rules: CP invariance implies the following selection rules for the photon- 
graviton conversion amplitudes [3]: 

• For a purely magnetic field e® couples only to e 1 - and e® only to £K 

• For a purely electric field e® couples only to E" and e® only to e 1 - . 
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This is borne out by the explicit calculation (see (26), (29)). 



Pair creation thresholds: In the purely magnetic case the amplitudes are real for small a>, 
since the magnetic field is not capable of pair production. The pair creation thresholds 
uo CT turn out to be identical with the ones for the corresponding photon-photon cases: 

2^1 + B, 

1 + y/l + 2B, (33) 
2. 



cr.scal 



cr,spm 
cr.spin 



to, 



cr,scal 



Calculable cases: The magnetic case is also much more amenable to an explicit 
calculation of the parameter integrals. In [TD] we have given a detailed analysis of 
the following regions in parameter space (with E = 0): 

• For photon/graviton energies below threshold the parameter integrals are suitable 
for a straightforward numerical evaluation. 

• For arbitrary u but small B the two-parameter integrals can be reduced to one- 
parameter integrals over Airy functions. 

• For u < uj cv and large B one finds the asymptotic behaviour 



(34) 



These leading asymptotic terms can be directly related to the corresponding UV 
counterterms, which is another property known from the photon-photon case [4"0] . 

In the zero energy limit, the amplitudes relate to the magnetic Euler-Heisenberg 
Lagrangians Cf^ spin (B): 

n®al,spin(^ — 0) B) — - ~^4~(^^ + ~0B^ ^™' spin ^ ' 
rr®|| t " r>\ 47TQ 1 d EH ' 

n sca l, S pin( u; — 0) B) — ^A^>^h ^scal,spin(-°) • 

III' JD OJD 

(35) 



The identities (35) have also been derived by Gies and Shaisultanov using a different 
approach [4*Tj . 
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5. Conclusions 



The calculation presented here is the first calculation of the photon-graviton vacuum 
polarization in a constant electromagnetic field, and also the first state-of-the-art 
application of the "string-inspired" worldline formalism to an amplitude involving 
gravitons. Although it was not possible here to go into detail, it should be emphasized 
that in this formalism this calculation is only moderately more difficult than the photon- 
photon polarization in the field. Moreover, we have also shown that the properties of 
the photon-graviton polarization tensor are very similar to the ones of the photon- 
photon one. We expect that even the graviton-graviton case will be quite feasible in 
this formalism. In a future sequel, we intend to analyze this case at the same level of 
the photon-graviton one, and to study the complete one-loop photon-graviton dispersion 
relations Q. 
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